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ABSTRACT 

CK Vul is classified as, amongst others, the slowest known nova, a hibernating nova, or 
a very late thermal pulse object. Following its eruption in AD 1670, the star remained 
visible for 2 years. A 15-arcsec nebula was discovered in the 1980's, but the star 
itself has not been detected since the eruption. We here present radio images which 
reveal an 0.1-arcsec radio source with a flux of 1.5 mJy at 5 GHz. Deep Ha images 
show a bipolar nebula with a longest extension of 70 arcsec, with the previously 
known compact nebula at its waist. The emission- line ratios show that the gas is 
shock-ionized, at velocities > 100 km s -1 . Dust emission yields an envelope mass of 
~ 5 x 10~ 2 M Q . Echelle spectra indicate outflow velocities up to 360 kms -1 . From a 
comparison of images obtained in 1991 and 2004 we find evidence for expansion of 
the nebula, consistent with an origin in the 1670 explosion; the measured expansion is 
centred on the radio source. No optical or infrared counterpart is found at the position 
of the radio source. The radio emission is interpreted as thermal free-free emission from 
gas with T c ~ 10 4 K. The radio source may be due to a remnant circumbinary disk, 
similar to those seen in some binary post-AGB stars. We discuss possible classifications 
of this unique outburst, including that of a sub-Chandrasekhar mass supernova, a nova 
eruption on a cool, low-mass white dwarf, or a thermal pulse induced by accretion from 
a circumbinary disk. 

Key words: Stars: CK Vul Stars: mass loss; Stars: evolution; Stars: binary; Stars: 
Post-AGB; Planetary nebulae 



1 INTRODUCTION 



The enigmatic CK Vul l|Shara fc Moffat! 1 19821 ) is cur- 
rently the oldest catalogued nova variables Its classifica- 
tion remains controversial. The star was first observed by 
Pere Dom Anthelme and was discovered independently by 



1 Of the eruptive variables, it is predated only by the well-known 
Galactic supernovae and by the LBV eruptions of P Cyg in 1600 
and 1655 



Hevelius in June/ July 1670. The observed eruption lasted 
for two years, during 1670-1672. It was discovered at 3rd 
magnitude, and faded over 100 days to below visual limits 
(magnitude 6). The following year it was recovered during 
rebrightening, reaching magnitude 2.6, 300 days after the 
first maximum. It again faded on the same time scale of 100 
days to below observational limits. It showed a faint tertiary 
maximum 600 day s after the first peak, at magnitude 5.5 
l|Shara et al.lll985 V The star has not been recovered since. 
Its light curve is unlike any catalogued novae, raising doubts 
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on its status. Interpretations include a hibernating nova, or 
a helium flash on a new white dwarf. 

The historic observations provided coordinates for CK 
Vul ac curate to within a few arcminutes. Based on the old 
charts, IShara fc Moffat] |l982t ) were successful in searching 
for the remnant. They discovered several nebulosities on 
their Ha+[N n] image, the spectra of which present reddened 
emissions of [Nil], Hn, [On] and [S n] . On their spectra, the 
flux of [N n] A 6584A is three times more intense than the 
Ha line, indicating an evolved, and perhaps hydrogen de- 
ficient chemical composition. No obvious exciting star was 
found in the field of CK Vul. A central star candidate was 
suggested bv lShara et al.l (Il985l), b ut shown to be misidenti- 
fied as such bv lNavlor et al.l (|l992| ); the latter authors were 
unable to identify, among the stars visible on their R im- 
age, an alternative central star candidate which could be 
responsible for the ionization of the observed nebulosity. 

The lack of bright candidates for the central star would 
constrain the possible companion of the white dwarf in cata- 
clysmic systems, assuming that CK Vul is an ordinary nova. 
However, the list of features uncommon for classical novae is 
extensive. The light curve during the outburst of Nova Vul 
1670 is unprecedented, showing variations which resemble 
a set of declines and subsequent recoveries. The rebright- 
enings may indicate constant intrinsic brightness with in- 
termittent dust formation, or a light echo. In the first in- 
stance, the object remained at high luminosity for over a 
year, with the fading(s) caused by intermittent extinction. 
In the second CELS6, cl dense medium needs to be present 
within a parsec of the source. The properties of the shell are 
quite unusual: a low outflow velocity of the ejected material, 
low density and hi gh ionized mass (a lthough the distance 
is very uncertain). IShara et al.l (|l985l ) discuss (and reject) 
some alternat i ve hy potheses concerning the nature of CK 
Vul. iHarrisonl l| 19961 ) suggests that CK Vul may be a Very 
Late Thermal Pulse event (VLTP) rather than a nova: this 
is based on similarities of the light curve to those seen in 
V605 Aql and V4334 Sgr (Sakurai's Object) (|Haiduk et all 
120051 ). Ie vans et all (|2002h have added more arguments for 
this reborn scenario for CK Vul. If the VLTP scenario for 
CK Vul is correct, it would be o nly the third object ob- 
served during a VLTP. lKatol (120031 ) suggested CK Vul to be 
a stellar merger. However, none of the presented scenarios is 
fully convincing: the observational data are limited, but the 
object seems to be unique. 

In order to test the VLTP hypothesis we searched for 
evidence of ionized material, since these stars reach high 
temperatures soon after the eruption and ionize a compact 
shell. In contrast, radio emission from classical novae fades 
within a few years after the eruption. To this end we con- 
ducted radio observations and also obtained deep images, to 
search for the ejecta, and for evidence of a fossil planetary 
nebula. 



2 OBSERVATIONS 
2.1 Radio emission 

The first radio observations of CK Vul were carried out with 
the VLA array on the 4 th of April 2005. The VLA was in the 
so-called 'B' configuration, with a largest baseline of about 
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Figure 1. The profile of the radio source, based on the higher 
resolution 8 GHz 2006 observations. 



10 km. The observations were done at a frequency of 5 GHz. 
CK Vul was observed for approximately 2 hours, divided 
into ~14-min scans, interspersed with observations of the 
phase calibrator (1922+155). A second set of observations 
was obtained on April 9, 2006, with the VLA in the 'A' 
configuration. After each ~10 min scan of CK Vul, at 5 GHz 
and 8 GHz, alternately, the phase calibrator 1922+155 was 
observed. 

1331+305 was used as the primary flux calibrator for 
both the 2005 and 2006 observations. Observational param- 
eters are summarized in Table [T] 

The data were reduced using standard procedures im- 
plemented in the AIPS package, by Fourier transforming the 
visibility data followed by the removal of the dirty beam pat- 
tern. The data were convolved with a synthesized gaussian 
point spread function (beam) of the same FWHM as the 
dirty beam: the beam sizes are listed in Table [1] Low spatial 
resolution images were also created, in order to search for 
confusing sources: only a few such confusing sources were 
found at 5 GHz and 8 GHz. Both uniformly and naturally 
weighted images were made, using the CLEAN algorithm. 
The RMS noise on the maps of CK Vul followed closely the 
theoretical thermal noise. 

Our 2005 observation reveals radio emission of 1.46 ± 
0.05 mjy (where the uncertainty is taken as the 3-<r noise on 
the map). The emission is consistent with a point source: the 
upper limit of the diameter, deconvolved from the beam, is 
0.5 arcsec. But the higher resolution 2006 observations find 
the object to be resolved. The profile of the source, and that 
of the beam are shown in Fig[T] The FWHM of the beam is 
0.26 x 0.23 arcsec, the convolved source 0.29 x 0.26 arcsec, 
and the deconvolved FWHM (gaussian fit) is 0.12 ± 0.01 x 
0.11 ± 0.01 arcsec. The fitting rules out a point source. 

The fluxes are listed in Table [T] they were measured 
both from gaussian fitting and from summing. The slight 
difference between the two epochs may be due to calibra- 
tion uncertainties or intrinsic variability. The peak inten- 
sity of the 8 GHz emission (natural weighting) is 1.25 ± 
0.02 mjy/beam, lower than the integrated flux. This also 
confirms that the source is resolved. The source is not cir- 
cularly polarized, but due to the faintness of the source the 
detection limit for polarization is only about 20 per cent. 
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Table 1. Observational parameters for the VLA observations (naturally weighted images). The radio source is at position (J2000) 
19 h 47 m 3 8 .o74 s + 27° 18' 45.16" 



Date Freq. t int Beam pa flux a 1922+155 flux 1331+305 flux 

FWHM deg mjy mjy/beam mjy Jy 



4 Apr 2005 5 GHz 2h 1.52"xl.36" -54 1.46 0.016 682 4.66 

9 Apr 2006 5 GHz lh 0.46"x0.38" +5 1.27 0.030 660 4.66 

9 Apr 2006 8GHz lh 0.26"x0.23" -1 1.53 0.023 651 7.34 



2.2 Ha observations 

2.2.1 INT 

The field of CK Vul was observed on two different occasions 
at the Isaac Newton Telescope (INT), located on La Palma. 
The first observations were carried out during a service night 
on Aug 4, 2004, and the second were obtained as part of 
the IPHAS survey, on July 12, 2005. IPHAS is a CCD Ha 
imag ing survey of the Northern Galactic plane l|Drew et al.l 
2005). This survey uses 120-sec exposures through the Ha 
filter, 30 sec through the Sloan r' and 10 sec through Sloan i 1 
(every sky position is covered twice). The deeper service ob- 
servations used three 600-sec Ha exposures, supplemented 
by three 90-sec exposures with the Sloan r' filter. All obser- 
vations used the wide field camera, which has a 30 arcmin 
field of view, mapped onto a mosaic of 4 CCDs. The pixel 
size is 0.3 arcsec. 

The observations were taken using the IPHAS observ- 
ing procedures and were reduced using the IPHAS pipeline: 
the procedure includes bias subtraction, flat-fielding, non- 
linearity correction and astrometric fit. The emission asso- 
ciated with CK Vul was centred on one of the four CCDs. 
The pipeline also obtains photometry of all point sources in 
the field, stored in a catalogue. The 2004 observations were 
taken under non-photometric conditions, with intermittent 
high cloud. 

The r' filter is used for the continuum subtraction, but 
it should be noted that its broad coverage (5500-7000A) also 
covers Ha. The Ha filter has a width of 95A, and includes 
the [Nil] lines. An emission line image was obtained from 
the difference of the Ha and r' images. 

The Ha and r' observations have different values for the 
image quality. A direct subtraction of the r' image therefore 
leaves significant stellar residuals-we note that the field is 
very crowded. Both positive and negative residuals are also 
caused by the fact that stars show variable r'-Ha colours 
due to the stellar Ha absorption, and by TiO bands for M 
stars. This can affect the r' — Ha colour by as much as a mag- 
nitude. To circumvent these problems, continuum subtrac- 
tion was done in a different way. We used the DAOPHOT 
program to obtain a list of stellar positions for the r' image. 
This list was subsequently used to find the counterparts of 
these sources in the Ha image. A point-spread function fit- 
ting was done to obtain the magnitude (in Ha), and to sub- 
tract these stellar sources. The same fitting and subtraction 
procedure was also carried out in the r' image. Finally, the 
residual r'-band image was subtracted from the residual Ha 
image. This procedure greatly reduces the stellar residuals 
but may subtract true Ha emission stars. We checked for 



the presence of such stars in the area surrounding CK Vul 
but did not find any. 

The astrometry of the image was re-calibrated using the 
positions of the stars in the Second US Naval Observator y 
CCD Astrograph Catalog (UCAC2 IZacharias et ai]|2004 ). 
The accuracy of the astromery, based on the RMS of the fit, 
is 0.1 arcsec. 

Below we refer to the final image as 'the Ha image' 
(adopting the name of the filter). However, the reader should 
be aware that the [N n] lines may contribute or even domi- 
nate the line emission within the filter pass band. 

2.2.2 WHT 

Ha images have previously been presented bv lNavlor et al.l 
(1992). The original data were kindly provided to us and 
re-reduced using IRAF. The observations were carried out 
at the William Herschel Telescope, on La Palma, on August 
10, 1991. They made use of a small CCD at the auxiliary 
Cassegrain port; the scale is 0.1 arcsec/pixel. The obser- 
vations in Ha were de-biased and flat-field corrected. The 
image was calibrated with respect to the positions of the 
stars on the INT image: the RMS of the relative fit is 0.05 
arcsec. 

2.3 Long-slit spectrometry 

The long-slit observations were obtained with the Ma nches- 
ter echelle spectrometer (MES: iMeaburn et al J 12003) com- 
bined with the f/7.9 focus of the 2.1 m San Pedro Martir, 
UNAM telescope, on 2006 July 12. This echelle spectrom- 
eter has no cross-disperser. For the present observations a 
filter of 90A bandwidth was used to isolate the 87 th order 
containing the Ha and [N n] emission lines. A 1024 x 1024 
SITE CCD with 24/im pixels was used. Seeing conditions 
were variable, with a mean of 2 arcsec during the observa- 
tions. 

Four slit positions across CK Vul were obtained. The 
first two positions employed an on-chip binning factor 2 x 
2 and a slit width of 150/xm, equivalent to llkms -1 and 
0.6 arcsec increments along the slit. However due to the 
faintness of the emission and the relatively poor seeing, the 
third and fourth positions were obtained with a 4 x 4 binning 
and a 300/xm wide slit corresponding to 20 kms -1 resolution 
and 1.2 arcsec increments along the slit for better sensitivity. 
The slit orientations of the latter are shown in Fig. [7] and 
the spectra in Fig. [8] The slits were oriented in all cases N- 
S and the integration times were 1800 seconds. The spectra 
were reduced and calibrated in wavelength against a Th-Ar 
lamp using standard procedures. 
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Figure 2. The INT Ho image of CK Vul, before continuum sub- 
traction. A logarithmic intensity scale is used to bring out the 
faint emission. The vertical bar is 30 arsec in length; North is at 
the top and East is left. 



2.4 Other data 

An earlier radio c ontinuum ob s ervat ion of CK Vul at 20 cm 
was performed bv lBode et aD (|l987h . on January 28, 1984. 
They found an upper limit of 1.52 mjy, with a spatial res- 
olution of 14 arcsec. Due to the low resolution, this limit 
includes both the compact source and the extended nebu- 
losity. 

Propo sed R-band and JHK counterparts are sum- 
marised in Evans et al] <|2002h . We consider these to be un- 
certain, referring to unrelated (field) stars. 2MASS finds all 
these to have the colours of normal stars. Observations in 
other wavelength regimes include flux measurements at 850 
and 450/^m and flux densities derived from re-analysis of the 
IRAS data (|Evans et alj|2002ft . Upper limits were found for 
MSX (8.3 pm) and IRAS 100 fim. 

GLIMPSE/Spitzer images of the field confirm that all 
2MASS sources show stellar Ray leigh- Jeans colours. This 
includes the star insid e the cavity sou th of the radio source, 
which was studied by iHarrisonl l)l996l l. 



3 RESULTS OF THE OBSERVATIONS 
3.1 Optical lobes 

The Ha image, derived from the Aug. 2004 (long) integra- 
tions, is shown in Fig. [2] It shows the dense star field , and 
the compact structure first found bv lShara et all l|l985h . The 
nebulosity appears as a partial ring, extending 15-arcsec 
across E-W. The star at the centre of the ring, originally 



suspected to be CK Vul, is now thought to be a field star. 
Faint nebular emission is lost in the dense star field. 

The corresponding continuum-subtracted optical image 
is presented in Fig. [3] The image shows that the previously 
known nebulosity is located at the waist of a much larger, 
bipolar structure. The bipolar lobes are edge-brightened: the 
filled emission is just visible in the northern lobe but is be- 
low the detection limit in the southern lobe. A number of 
small condensations are visible, two of which delineate the 
tips of the lobes. A bright condensation is seen on the SE 
lobe. The angular distance between the tips of the lobes 
is 70 arcsec. The filter transmission includes both Ha and 
[N n] , and either line may contribute in different areas of the 
structures. 

The central region is shown in more detail in Fig. [3)3. 
A funnel-like structure is seen, with one side much brighter 
than the other. (But whether the structure is in fact a cavity 
is not known.) Four other knots are visible: e ach of these 
compo nents are numbered on the image as in IShara et al] 
l|l985h . An arc connects knots 2 and 4, but it is not part of 
a closed ring. The three-dimensional structure cannot easily 
be deduced f rom t he image. In addition to the numbering of 
IShara et al] |l985|), we number (6) the central diffuse source 
and (7) the more diffuse emission between knots 4 and 5. 

We note that a stellar source on the western ring was 
subtracted in the continuum removal process. This left the 
blank residual below knot 3. Some compact Ha emission was 
likely subtracted in this process, as suggested by the residu- 
als. This problem occurs when a star is much brighter than 
coincident Ha emission, but dos not affect other locations 
within the nebulosity. Fig. |4j\ shows the inner region of r'- 
band INT image, where this star is visible to the West of 
the radio source. 



3.2 Radio source 

The radio emission is shown in Fig. |4j\ superposed on the 
INT r'-band image. The brightest emission-line regions are 
faintly visible in the r'-band, because the wavelengths of 
the emission lines fall within the filter transmission curve. 
The radio source is located between knots 1 and 6 (Fig. 
[3|. Three stars located in the region below the funnel could 
be candidates for the central star, but they do not coincide 
with the compact radio source. The radio source also does 
not coincide with any of the Ha+[Nn] bright knots. 

An accidental superposition of a background radio 
source is improbable, due to the excellent coincidence of the 
position of the radio emission with the nebula. Moreover, 
only one other source with comparable flux to the one re- 
ferred to CK Vul was found in the full primary beam field 
(FWHM) of 9.0 arcmin. 

The radio source is close to the brightest part of the 
elongated knot T'. Fig. [4)3 shows the WHT image, which 
has better spatial resolution. A faint elongation or separate 
source is indicated just below knot T; on the INT image this 
merges with the knot because of the worse seeing. The faint 
elongation is within 0.5 arcsec of the radio source (right-most 
plus sign in Fig. [4)3). It has a very faint r'-band counterpart, 
confirming it is an emission line region. To classify this ex- 
tension as a separate source would require either higher an- 
gular resolution, or kinematic data. Given the extent of the 
emission line regions around the radio source, the likelihood 
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Figure 3. The continuum-subtracted INT Ha image of CK Vul. 
North is top, east is left. A (Top): Image encompassing all de- 
tected line emission. The contrast is set to bring out the faint 
bipolar lobes. The vertical bar is 30 arsec in length. B (Bottom): 
The central region , corresponding to the emission detected by 
Sh ara et all l|l985h . The linear contours and the grey scale refer 
to the same image. The peaks on either side of the apparent 'fun- 
nel' are 5 arcsec apart. The numbe rs 1 to 5 correspond to the 
knots detected by IShara et al.l lll985h . Numbers 6 & 7 have been 
added in this paper; see text for discussion. 



of a chance superposition is rather high. However, it is our 
only -but speculative- candidate for an optical counterpart. 

The radio source is situated in a cavity close to the 
brightest part of the nebulosity. Nebulosity located further 
from the radio source shows a tendency of decreasing bright- 
ness with distance from the radio source. This indicates that 
the radio source is associated with the current energy source 




Figure 4. The radio source superposed on the optical images. 
North is top, east is left. The position of the radio source is 
19 h 47 m 38.074 s +27° 18' 45.16". The knots are indicated by their 
numbers. A (top): superposed on the r'-band INT image. The im- 
age is 21.5 arcsec across horizontally. Radio contours are at 0.075, 
0.15,0.3,0.6,0.9,1.2 mjy/beam, from the 2005 VLA observations 
(natural weighting). The oval indicates the size of the beam; the 
black bar is 5 arcsec long. B (bottom): superposed on the WHT 
Ha image. The right-most plus sign indicates the position of the 
radio source, the left-most plus sign gives the geometric centre of 
the bipolar lobes. 

of the nebula, and reveals the position, if not the identifica- 
tion, of the central object. 

The radio source is to a very good (1-arcsec) accuracy 
located at the centre of symmetry of the tips of the bipolar 
nebula. This point is indicated by the left-most plus sign in 
Fig. 2)3. This also supports the identification of the radio 
source as the central object. 

The INT data show that at the position of the radio 
source, there is no optical counterpart down to a conserva- 
tive limit of r' > 23. The i-band IPHAS image is less deep, 
and gives a limit of il > 20. 

3.3 Expansion 

The WHT image compared to the INT image allows us to 
search for temporal variations in the nebula. The WHT data 
were taken in 1991 and the INT data in 2004. The 13-year 
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Figure 5. The WHT Ha from 1991 (grey scale) together with 
the 2004 INT Ha image (contours at 3 and 10 a). Two compact 
components show significant proper motion. The arrows link the 
radio sources to these components: the direction is consistent with 
the observed proper motion. 



separation represents 4 per cent of the elapsed time since 
the AD 1670 eruption. 

In Fig. [5] the grey scale shows the WHT image, and the 
INT Ha (including continuum) image is overlaid as contours. 
The WHT image does not cover the full region of the nebula 
(which was not known at the time). The images are aligned 
on the stellar positions. The relative (rms) accuracy is 0.05 
arcsec. 

Several emission-line nebulosities have changed posi- 
tion, in particular knots 4 and 5. (The object to the West of 
knot 3 has moved north: this is interpreted as a high proper 
motion star.) The change in position is away from the central 
radio source, as indicated by the arrows in Fig. [S] The posi- 
tions of knots 4 and 5 are given in Table[2] The indicated rel- 
ative motion from 1991 to 2004 is (/x a ,/Lt^) = (-0.16, -0.46) 
arcsec, and (n a ,jj,s) = (0.36,-0.67) arcsec, for knot 4 and 
5, respectively. 

The original positions of the two components are ob- 
tained by extrapolation back to 1670. The extrapolated po- 
sition of knot 4 is coincident with the position of the radio 
source, assumed to represent the central star. The agreement 
is worse for knot 5, possibly due to its extended, irregular 
structure, which is better resolved in the WHT image than 
in the INT image (Table Residual image distortions may 
also be present in the WHT images. 

The observed motion of these two components, located 
outside the inner region of the nebula, implies that the ob- 
served lobes are indeed the remnants of the 1670 event. It 



Table 2. Centroid positions (J2000) of the two components in- 
dicated in Fig. \5\ The accuracy of the determination of centroids 
is typically about 0.02 arcsec, but for the component 5 at 2004 is 
a few times worse. 
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Figure 6. The residuals of the scaled IPHAS image and the WHT 
Ha images, versus the position of the assumed expansion centre. 
The position of radio source is at 0. 



also confirms our interpretation of the radio source as the 
eruptive object. 

To improve the sensitivity to expansion, we subtracted 
the IPHAS image from the WHT observations. The sub- 
tracted image showed residuals at the position of the two 
outer, moving knots, the central brightest nebulosities, the 
faint arc and other nebular structures, present on the Ha 
image. Applying a constant scaling factor prior to the sub- 
traction, to reduce the image scale of the IPHAS image, 
greatly reduced these residuals. 

Figure [Tj] shows the residuals of the two knots, and the 
brightest part of the inner nebula, for different positions of 
the centre of the expansion. The best fit was obtained for 
the assumed expansion centre lying within 2 arcsec of the 
radio emission. The constant scaling factor assumes that the 
material seen in the field of CK Vul is coming from the 1670 
eruption event and is flowing ballistically. 

Unfortunately, the tips of the lobes of the outer nebula 
are outside of the WHT frame. However, the tip of the south- 
ern lo be is just visible on the plate published in lShara et alj 
(1985) (their Fig. 8). Comparison shows that this tip has 
moved by roughly 2.5 arcsec in the intervening 20 years. 
Within the uncertainties, this proper motion is consistent 
with a 1670 ejection event. 
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Table 3. Line fluxes for the different knots, as measured from the 
INT image. For knot 3, the flux includes a non-subtracted star. 
The values for knots 6 and 7 are also uncertain. The uncertainty 
in the flux conversion factor is estimated at 10 per cent. The [N o] 
6584/Ha ratios are indicative only. 



Knot 


Flux 

[10~ 15 erg cm -2 s" 1 ] 


v hcl 

[kms ] 


[Nil] /Ho 


1 


43 






2 


28 


-80 


1.5 


3 




+0 


1 


4 


4.7 


+60 


> 10 


5 


3.6 






6 


1.7 






7 


2.6 






southern tip 




-175 


> 2 



3. 3. 1 Brightness 

We measured fluxes for the different knots as follows: The r' 
magnitudes of 12 field stars were obtained from the IPHAS 
catalogues. The magnitudes were converted to Bessell niters 
on Landolt standard stars using the colour transformations 

T - 7? L an = +0.275(7? - J) L an + 0.008 

(r - i') = 1.052(7? - 7) Lan + 0.004 

Th e magnitudes were co nverted into flux units per A using 
the iBessell et ail (|l987f l coefficients for R. The integrated 
in-band flux was obtained by multiplying by the equivalent 
width of the filter (80 A); this was divided by the data count 
(ADU) for each star to get a flux per data count. We find 
a conversion value of (94 ± 6) x 10~ 20 erg cm" 2 s" 1 per data 
count. This factor was applied to the data counts for each 
of the knots in the continuum-subtracted image. 

The results are listed in Table [3] The line fluxes contain 
contributions from Ha, and from the two [N n] lines: the 
filter transmission is within 5% of the peak transmission 
for each of the lines. Based on the spectroscopic data, Ha 
accounts for typically half to a quarter of the measured line 
flux, but for knot 4 we do not detect Ha. 

The changes between the two INT images, measured in 
the same way, indicate a 0.08 mag (approx. 10%) change in 
the brightness of the knots, but this varies between knots by 
more than 10%. The standard deviation for the background 
stars alone is 0.05 mag, so that we do not see conclusive 
evidence for brightness evolution. 

The brightest knots are those closest to the radio source. 
This suggest that their excitation arises from a central 
source. 

3.4 Spectroscopy 

Our echelle spectra show [Nil] 6584/Ha ratios between 1 
and 5. Individual values are listed in Table [3] Because of 
the low S/N, these need confirmation. Knot 4 has no Ha 
detected and could have a very high ratio. 

These ratios can be compared t o other values reported 
in the literature. IShara et al.l l|l985h analyze spectra of the 
three brightest nebulosities, numbered 1-3. They find con- 
sistent line intensities, with remarkably strong [S n] emission 




Figure 7. Positions of the two long slit echelle spectra discussed 
in the text. The scale bar is 20 arcsec long. 



comparable to Ha, and a ratio [Nn]6584/Ha « 3. The [Sn] 
6716/6730 line ratio is approximately two. [O ill] is also seen, 
several times brighter than H/3, but weakened by the extinc- 
tion. 

A spectrum of a faint knot is reported bv lNavlor et al.l 

ljl992h . The knot has [Nil] 6584/Ha « 1, and a [Sn] 
6716/6730 ratio of unity. It is not clear which knot this rep- 
resents: their slit orientation at position angle 45 degrees 
covered knot 4 and 6. Our much higher line ratio for knot 4 
suggests Naylor et al. observed knot 6. 

A spectrum taken with a 4-arcsec slit ori e nted perpen- 
dicular to this direction is reported bv lCohenl (|l985l ). Their 
spectrum shows intermediate values for both the N and S 
line ratios compared to the previous two papers. 

3.5 Kinematics 

The two long-slit spectra with better signal to noise are 
shown in Fig. [S] only signal from the relatively bright knots 
in the core was recorded. The top and middle panels in Fig. 
[8] show both the Ha and [Nil] 6584A emission lines for each 
slit position. The locations of the slits are shown in Fig. [7] 
The first slit (top panel) intersects knot 2 and a section of 
the southern arc, whereas the second slit (middle panel) lo- 
cated to the West of the previous one, intersects knot 3 and 
the edge of knot 4. The lower panel shows the correspond- 
ing P-V arrays for the [Nil] 6584A line with heliocentric 
velocity scales. 

The spectra show that the [Nil] 6584A line is consis- 
tently stronger than Ha. Although Ha is barely detected, 
the S/N is enough to show that its line profile follows the 
same wide structure of [N n] 6584A. The line profiles are 
peculiar in the sense that the emission seems to arise only 
from localised regions and with a wide velocity range, of 
the order of 350 kms" 1 , somehow akin to tho se observed in 
supernovae blasts (e.g. lRiesgo fc L6pezll20o"5h . 

The brightest region of each line profile yields the fol- 
lowing heliocentric velocities: for knot 2 we find an approx- 
imate velocity of —80 km s" 1 ; for knot 3, km s" 1 ; for 
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Figure 8. Two long slit echelle spectra. The upper frames show 
the full order, covering Ha and [N II] 6584A. The two lower frames 
show only the [N II] line. Lower left frame: knot 2 (upper) and the 
southern arc (lower emission). Lower right frame: knot 3 (upper) 
and knot 4 (lower emission). 



knot 4, +60 km and for the southern arc, +70 km . 
From these line profiles the systemic velocity cannot be de- 
termined with confidence. However, since the observed he- 
liocentric velocities are heavily blue-shifted for knots 2 & 
3, it is reasonable to assume that there we are detecting 
material that is mainly coming towards us whereas the op- 
posite is observed in knot 4 and the southern arc. From this 
kinematic behaviour and the morphology of the region (see 
Fig. [3} a possible interpretation is that of a rapidly expand- 
ing (and disrupted) torus or equatorial structure, centred at 
roughly —50 km s _1 and where knots 2 & 3 are part of the 
foreground (approaching) equatorial structure and knot 4 
and the southern arc are part of the background (receding) 
part of the torus. 

The tip of the southern lobe was covered by two of the 
slit settings. There is a faint trace in both spectra of a ve- 
locity feature at the location of the tip with a velocity of 
about —175 km s" 1 . The sign and magnitude of this velocity 
is consistent with the bipolar lobes being oriented perpen- 
dicular to the equatorial structure discussed above and with 
the northern lobe moving away and the southern one coming 
towards us. 

The kinematics and structure of the line profiles dis- 
cussed here indicate that shocks are a main, if not the only, 
contributor to the excitation of these emission lines, in agree- 
ments with the arguments of the preceding section. 



IShara et"afl (|l985l ) locate the nebula bey ond an extinc- 
tion l ayer, located at 550 ± 150 pc distance. IWeight et al.l 
(|l993t ) deduce an upper limit of 2kpc: they detect two CO 
clumps, with kinematic distances of 500 pc and 2 kpc. The 
two clumps together would result in E(B — V) = 1.25, sug- 
gesting that CK Vul must lie in front of the more distant 
clump. 

We have assessed the IPHAS photometry of point 
sources in the vicinity of CK Vul to constrain the distance- 
extinction relation in this field. Within a field of 4.5' x 4.5', 
we find a few essentially unreddened main sequence stars, 
and then a dearth of stars to A v ~ 1.2 (E(B - V) ~ 0.4). 
This jump in extinc tion may rep r esent the dust layer at 
550 pc described by IShara et al.l (|l985l ). We have cross- 
identified a group of 10 early A stars (presumed to be 
dwarfs) for which there are 2MASS JHK measurements as 
well as IPHAS r' , i' and Hcv data, and used them to obtain 
a sample of reddenings and distances consistent with their 
likely spectral types. These early A stars range in bright- 
ness between r' ~ 14 and r' ~ 18, and are consistent with 
extinctions in the range 0.9 ± 0.2 < Ay < 5.4 ± 0.5. The 
deduced distances fall in the range 2-4 kpc. The reddening 
does not increase in any organised manner with distance, 
and the extinction layer(s) appear to be patchy. 

We will here use a distance of 600 pc. The extinctions 
indicate an upper limit for the distance of CK Vul of 4 kpc. 



4.2 Velocities 

The radial velocity range is about 350 km s -1 . The highest 
blue-shifted velocities are seen just to the north of knot 2 and 
may represent the inner edge of the bipolar lobe, oriented 
close to the line of sight along the edge of the torus. We do 
not see a similar extended velocity structure at red-shifted 
velocities. This may be due to the placement of the slits, 
or may indicate some internal extinction at that location. 
The southern tip is blue-shifted at — 175 kms -1 : similar to 
the extreme velocity near the torus, but likely less inclined 
towards us. 

The proper motion of the knots 4 and 5 are around 
37 mas/year and 59 mas/year respectively. For the outer 
tips of the bipolar lobes, the radial distance to the central 
radio source of 35 arcsec indicates a proper motion of 105 
mas/year. Our measured proper motion of the southern tip 
corresponds to 120 ± 25masyr _1 . For a distance of 600 pc, 
the respective tangential velocities are 106 and 169 km/s 
for the knots, and 340kms _1 for the tips of the bipolar 
lobes. Assuming a systemic velocity «h e i = — 50kms _1 , the 
space velocity of the tip of the bipolar lobe becomes w e j ~ 
360 kms -1 . 

The bipolar flow is inclined by 20 degrees with respect 
to the plane of the sky, from these numbers. 



4 PHYSICAL PARAMETERS 
4.1 Extinction distance 

The spectra of IShara et alj (1 19851 ) show a strong Balmer 
decrement, corresponding to E(B — V) = 0.82 ± 0.23, or 
A(V) — 2.5). Based on the weighted mean of this and the 
extinction towards a distant field star ( E(B — V) = 0.7 ± 
0.1), they adopt A(V) = 2.2 ± 0.3. 



4.3 Extended nebula 

The nebula is detected only in ionized lines and in dust. It 
is possible that emission lines reveal only part of the nebula, 
with the remainder not ionized or too faint. Parameters of 
the extended nebula m ust therefore depend on estimations. 

IShara et all l|l985l) d erive n e < 10 cm , based on the 
[S n] line ratio in knot 1. iNavlor et ail l|l992l) find a higher 
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density n c ~ 10 3 cm -3 , based on the [Sn] line ratio of a 
faint knot (possibly knot 4 or knot 6). If we assume that 
the system of knots 1,2,3,4,6 represents a partially ionized 
torus with density n = 10 2 cm -3 , diameter of 10 arcsec, 
thickness and height of 5 arcsec, the total mass becomes 
M ~ 2 x 10~ 5 Mq. If, on the other hand, we assume that the 
total extent is given by the bipolar nebula (35 arcsec radius), 
and assume a typical density of n — 10 2 cm -3 , the total 
mass becomes M ~ 2 x 10" 2 M Q . The dust mass derived 
bv lEvans et all (|2002h gives M~5x 10" 2 M Q for a typical 
dust-to-gas ratio of 100 (assuming hydrogen-rich gas). Thus, 
a model of a cloud of radius 30 arcsec and average density 
n = 10 2 cm" 3 is c onsist ent with observational constraints. 

IShara et all (|l985l ) suggests that the extended nebula 
was ionized during the 1670 eruption and is now recom- 
bining. This assumption of a slowly recombining nebula 
presents some problems. N + recombines much faster than 
hydrogen, so the high [Nll]/Ha ratio would require that no 
significant recombination has yet occurred even in the denser 
knots. For [Nil], at a density of 100 cm" 3 the reco mbination 
time scale is 75 yr (jLechner fc Kimeswen gcr 2004). Densities 
of n c <i 10 2 cm -3 require contin uous ionizati on. 

The line ratios listed by Sh ara et al ] l|l985l ) are m 
go od qualitative agreement wit h the radiative shock models 
of iMcKee fc Hollenbachl (|l980h . although Ha is somewhat 
fainter than seen in those models. In particular, the [S n] 
6716+6730/Ha ~ 1, [Nl] 5200/H/3 ~ 1 and the strength 
of the [Nil] 658 4A line are indicative of shock velocities 
v ~ 100 km s" 1 . (jShara et all (|l985l ) do not list or identify 
the [N i] 5200A line, but it is apparent in their spectrum.) 
The presence of [Oiii], if shock-ionized, indicates shock ve- 
locities above lOOkms - . 



5 THE COMPACT RADIO SOURCE 
5.1 Emission mechanism 

The fact that we have a compact but resolved source sug- 
gests thermal free-free emission for the origin of the ra- 
dio emission. This is supported by the observed brightness 
temperature. At 5 GHz, the observed diameter (0.11 arcsec 
FWHM of a deconvolved gaussian) gives a brightness tem- 
perature of T b = 6000 K. At 8 GHz we find T b = 2000 K. 
Assuming the electron temperature is T c — 10 4 K, the opti- 
cal depth at 5 GHz is r « 0.9. This predicts at 8 GHz that 
Tb ps 2700 K, in reasonable agreement with observations. 
The emission is optically thick at 1.4 GHz, expla ining the 
non-detection at this frequency l|Bode et al.lll987h . 

The lack of an infrared counterpart presents a problem, 
however. Optically thin free-free emission predicts approxi- 
mately 0.6 mjy at 8/xm, assuming a v~ 01 index. The upper 
limit for a counterpart to the radio source at 8/^m is around 
0.2 mjy. The only possible non-stellar source is a possible 
faint source 2 arcsec south-east of the radio source, which is 
absent from the Spitzer images at shorter wavelengths. The 
lack of an 8/im detection can be explained by very high ex- 
tinction. However, it also suggests to explore different emis- 
sion mechanism. 

Small grains can emit appreciable radio emission 
through their rotation. The ele ct ric d ipole emission is 
described in iDraine fc Lazarianl (| 19981 ) and references 



therein. Using emissivitie s per hydrogen atom tabulated in 
IDraine fc Lazarianl (1 19981 ) , we find that our radio flux would 
require ~ 3 x 10 -3 Mq of gas (assuming normal composition 
of the gas). This is a plausible value for the radio source, 
and would give densities around n ~ 10 7 cm~ 3 . However, 
this emission produces a positive spectral index between 5 
and 8 GHz, while we find a flat flux distribution. There is 
also no evidence for the required very small grains in the 

infrared spectral energy distribution. 

Synchrotron and cyclotron emission l|Linksvl 19961 ). from 
relativistic electrons moving in a magnetic field, give bright- 
ness temperatures much higher than we observe. However, 
if our radio source is found not to be resolved, this becomes 
a possibility. 

5.2 Density and mass of the radio core 

If the radio emission is due to free-free emission, we can 
derive some parameters of the originating nebula. Using the 
observed flux and assuming r = 0.9 and T c = 10 4 K at 
5 GHz, the emission measure becomes EM = n 2 / = 8.1 x 
10 7 pc cm -6 where n e is the electron density and I the depth 
of the emitting region. The latter is taken as the same as the 
FWHM of the radio source (assuming symmetry) , or 66 AU 
at 600 pc. 

We arrive at a density of iV e = 5 x 10 s cm -3 . The mass of 
the ionized region becomes Mi ~ 4 x 10 -7 Mq. This assumes 
a H-rich chemical composition: H-poor gas requires more 
mass per electron. 

For a filling factor equal to unity, a relation between 
the hydrogen line flux and the optically thin free-free Ssghz 
flux density is given bv lPottasc h (1984) assuming hydrogen- 
rich material. This gives an expected Fo(Hq) ~ 1.5 x 
10" 12 ergs -1 cm" 2 . The IPHAS Ha filter has an equiva- 
lent width of a pproximately 80A and a peak transmission 
of 90 per cent ijDrew et al.ll2005l ): the unreddened magni- 
tude of the central object in the Ha filter should have been 
m ~ 12.3. The (undetected) object is much weaker. The 
emission located 0.5 arcsec away from the radio source is 
1.6 x 10 _15 ergs _1 cm -2 , three orders fainter. This requires 
An ^ 7.5, corresponding to Av 1? 10 for a standard extinc- 
tion law and pure Ha emission. 

The [Nil] 6584/6548A lines are within the IPHAS Ha 
filter transmission. However, these lines have a critical den- 
sity of n cr = 9 x 10 4 cm -3 . We derive an electron density 
of n c = 5 x 10 5 cm~ 3 , suggesting that these lines do not 
contribute in the central region. 

The lack of a clear Ha counterpart to the radio source 
can also be interpreted as evidence that the core is H-poor. 
Hydrogen-poor gas tends to have lower electron tempera- 
tures, as the cooling per nucleus is more efficient. We ran 
some Cloudy models to obtain in dicative values . The mod- 
els are similar to those described in iHaiduk et al] (2005) , but 
use a PG1159 stellar model (|Rauchll2003l ). with luminosity 
L — 0.5 Lq, a density n e = 5-7 x 10 5 cm -3 (depending on 
stellar temperature T^ff ) and an inner radius of 13 AU. For 
T cft = 60, 80, 100, 120 kK, we find T c = 8100, 8400, 8700, 
9200K, respectively. These are lower than corresponding val- 
ues of H-rich gas, but not dramatically so. The physical 
cause of the relatively high temperatures is the suppression 
of the forbidden lines at these high densities, which reduces 
the cooling rate. The model temperatures are within the 
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range derived from the radio free-free spectrum. But in the 
absence of measured line ratios, the models are exploratory 
only. 

5.3 Luminosity 

The recombination timescale of the radio core is less than a 
year. The central region is therefore continuously ionized: al- 
though no central star is visible, its continuing existence can 
be inferred. Recombining gas can also be ruled out because 
it would show a very low electron temperature T e ~ 10 3 K, 
much lower than the observed brightness temperature. 

The luminosity of the ionizing star can be estimated as- 
suming that all photons above the Lyman limit ionize hydro- 
gen, and that the object emits as a black body. The number 
of ionizing photons reaches a maximum at a temperature of 
about T* = 7 x 10 4 K. For this temperature, the required 
black body luminosity is L* = 0.5 Lq. For other tempera- 
tures within the range 40-100 kK, the luminosity is within 
a factor of 1.5 of this value. Stellar atmosphere models re- 
quire a little higher luminosity as a black body emits more 
ionizing photons than a realistic atmosphere of the same 
temperature (e.g., at T = 50 kK, a Kurucz model gives 10 
% fewer H-ionizing photons than does a black body). Geo- 
metrical considerations also indicate that this luminosity is 
a lower limit, and the real luminosity could be a few times 
higher. 

Such a luminosity, combined with a high stellar temper- 
ature, gives an unreddened magnitude of order Vo ~ 17. No 
star is detected near this magnitude. The ROSAT all sky 
survey did no t detect any source at the position of CK Vul. 
ISevellil (|2004l ) reports extreme faintness of CK Vul in the 
SWP IUE spectra; the aperture of the instrument included 
the position of the radio source. The lack of detection may 
be due to a higher stellar temperature together with high 
obscuration. 

For comparison, at the adopted distance, the peak 
absolute magnitude of the AD 1670 eruption was around 
M (V)o ~ —9 to —8. The current luminosity implied by the 
radio source is a factor of 10 5 fainter. 



5.4 Dust 

The infrared spectrum of CK Vul l|Evans et alj r2002l shows 
two components, with dust temperatures around 550 K and 
25 K. Their two-component dust model yields approximate 
dust masses of 1.5 x 10 -9 and 5 x 10 -4 M . The mass of the 
first component agrees well with the derived ionized mass, 
assuming normal dust-to-gas ratios. The model was based on 
IRAS flux densities: the large IRAS apertures include both 
the core and the extended nebula, and the second component 
may be identified with the extended nebula. 

The new high-resolution Glimpse limit is a factor of 
100 below the IRAS 12/xm flux. We ran dust models with 
parameters n = 2 x 10 14 cm, r = 10 15 cm, density at inner 
edge 5 x 10 s cm -3 , density distribution n oc r , and a star 
with T = 6 x 10 4 K and L = 0.5 Lq. The (assumed H-rich) 
shell has a mass consistent with the ionized mass. For silicate 
dust, this yields predicted flux densities of F$ = 0.1 mjy and 
F 12 ~ 0.6 mjy, well within observationa l limits. Carbon-rich 
dust (as assumed bv lEvans et al.l ((2002) for their hotter dust 



component) would be several times brighter at 8 /xm. These 
exploratory models differ from lEvans et all (|2002h mainly in 
that they give lower dust temperatures (200-300 K). 

The peculiar light curve with several dips may be due 
to epochs of dust formation. In support of this, Hevelius 
comments that during t he second peak ( April 1671), the 
source appeared reddish (|Shara et al.lll985l ). The brightness 
of the initial eruption suggest that it did not suffer the line- 
of-sight extinction of the current radio source, and that the 
dust formed afterwards. 



6 THE NATURE OF CK VUL 
6.1 Observational constraints 

The current luminosity and high temperature places the cen- 
tral object high on the white dwarf cooling track, e.g., a 1 Lq 
white dwarf has T e g « 65 kK and a cooling age of t c ~ 10 6 yr 
(|Bloeckerll 19951 ). 

The compactness of the radio source implies it is not 
expanding, as the expansion velocity would be less than 
lkms -1 , or below the escape velocity. We can therefore as- 
sume that the radio nebula is part of a stable structure, such 
as a rotating disk, and could predate the eruption. The bipo- 
lar lobes are similar to those found in some young planetary 
nebulae and post-AGB stars. Such lobes can form where a 
spherical f ast wind meets an existing equatorial obstruction 
(|lckdll987l ). 

An area of low du st emission surrou nding CK Vul is 
seen in IRAS images l|Evans et al.ll2002l) . This resembles 
cases where the ISM has be en swept up by a stellar wind 
i|Ziilstra fc Weinbergerl l2002 ) and suggests CK Vul has gone 
through phases of high mass loss over the past ~ 10 5 -10 6 yr. 

Together, these constraints favour an evolved object, 
and the morphology favours a binary object. Two possible 
interpretations are a nova, and a thermal pulse. 

One should consider whether a supernova classification 
can truly be ruled out. The [Sii]/Hq ~ 1 ratio is diagnostic 
of supernova remnants. Also, the current luminosity is a 
typical value for a pulsar, a nd the radio spectr um is similar 
to pulsar wind nebulae (e.g. iTemim et al.ll2006i ). The object 
H2-12, a knot inside the Kepler supernov a remant, shows 
simil ar line ratios and velocities to CK Vul (|Riesgo fc LopeJ 
120051 ). However, there is no X-ray source at this position, nor 
is there significant radio emission from the extended nebula. 
The bipolar symmetry of CK Vul differs from the structures 
seen in supernova remnants. We conclude that a supernova 
classification is unlikely. 

A 'gentle' supernova (lToutll2006h can be considered. An 
accreting helium white dwarf may ignite its helium at masses 
far b elow the Chandrasekhar mass |Wooslev fc Weaver] 
1 19941) . The lack of observational evidence for sub-luminous 
supernovae suggests that the ignition under such conditions 
may be non-explosive. However, such ob jects may also have 
been classified wrongly. The models of Iwooslev fc Weaverl 
l|l994h predict short-lived maxima, however. 

The observatio ns do not cl early support the possibility 
of a stellar merger (|Katoll2003j ): a merger would be able to 
explain the large drop in luminosity since the eruption, but 
would lead to a bloated, cool star. The radio source requires 
an ionizing star. 
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6.2 A light nova 

Classical novae are dominated by massive white dwarfs: 
these show frequent eruptions, but consequently eject rel- 
atively little mass. A lower-mass white dwarf with a low 
accretion rate can eject the most massive shell. However, 
such eruptions are rare as a consequence of the much longer 
accretion time scales. The current luminosity of CK Vul cor- 



responds to an accretion rate of A/ a , 



10" 



'M yr _1 . The 



peak temperature, for a thermal-nuclear runaway on a low- 
mass white dwarf, may not even reach 10 s K limiting the 
nucleosynthesis wich occurs (|Starrfieldl2007h . 

A large set of mode ls covering a range of p arame- 
ter space is presented by iPrialnik fc Koveta (| 1995T) . Some 
of these approach the constraints set by CK Vul. For in- 
stance, their model with Mwd = 0.65 Mq, Twd = 3 x 10 7 K 
and M = 10 -10 Mq yr -1 yields an ejection velocity of 
200 km s" 1 , eruption amplitude of 16 mag, mass- loss time 
scale of 480 days (during which time the nova remains 
bright), and a very high peak luminosity of L — 8 x 10 4 Lq 
(Mbol = —7.5). The recurrence time scale of this model is 
10 6 yr. 

This particular model gives an ejecta mass of 10 -4 Mq. 
The highest p rediced ejecta mas s for a light nova is only 
7 x 10" 4 M Q l|Yaron et all 120051 ). On the one hand, this 
is three orders of magnitude higher than found in near- 
Chandrasekhar-mass novae. On the other hand, it is one to 
two orders of magnitude below the derived nebular mass of 
CK Vul (and below the mass in the dust alone). It is how- 
ever possible that part of this is swept-up material which 
predates the eruption. 



6.3 A thermal pulse 

Helium ignition on the white dwarf cooling track has been 
observed in V605 Aql (IClavton et al.ll2006l) and V4334 Sgr 
(Sakurai's Object : lAsplund et al.lll999l : lHaiduk et al.ll2005l : 
lEvans et aljEjbl) . Their light curves show similarities to 
CK Vul: a long lived maximum at M ~ —5, with a brief 
phase of very high, dusty mass loss. Three further objects 
are known to have experienced a VLTP event in their evo- 
lution, as shown by the presence of hydrogen-poo r cen- 
tral regions within a planetary nebula (|Ziilstrall2002l ). Out- 
flow Velocities are also si milar, with observed values are 
~ 25 ferns -1 for V605 A ql (|Pollacco et al.lll992T ) and V4334 
Sgr (|Kerber et al.ll2002f i. 

If CK Vul is a VLTP, it would fall on an evolution- 
ary sequence V4334 Sgr (Sakurai's Object: 10 yr), V605 
Aql (90 yr), IRAS 15154-5 258 (~ 10 3 yr ) and the two more 
evolved objects A30/A78 jZiilstral 12003 ). The first two are 
still hidden inside a dense cocoon of dusty ejecta, but like 
CK Vul the central object is seen in radio emission. There 
are some morphological resemblances with CK Vul, with the 
ejecta in these objects showing equatorial torii and polar 
knots. The origin of the morphologies are not known. 

There are however also notable differences. The radio 
source in CK Vul is much more compact than expected: 
in V605 Aql, which is at a larger distance, the 0.3-arcsec 
dense shell is expanding and there is no evidence for a non- 
expanding component. The ejecta in both VLTP objects are 
almost completely H-poor. In contrast, most knots in CK 
Vul show significant Ha components, even if fainter than 



expected. (We do not detect Ha emission in knot 4, which is 
likely H-poor. It is possible that the other knots have swept 
up ambient hydrogen.) CK Vul has a much lower current lu- 
minosity than the other five objects. Finally, fossil planetary 
nebulae are seen in all five objects, but CK Vul only shows a 
surrounding hole in the ISM as evidence for past mass loss: 
this suggests a much longer post-AGB time scale before its 
eruption. 

VLTP events may show a characteristic double loop in 
the HR diagram, ov er a time scale of a few hundred years 
jHaiduk et al.ll2005l ). Direct observational evidence for a sec- 
ond loop i s still lacking and the VLTP sc enario is not well un- 
derstood ( Mill er~"Bertolami et al.ll2006l). Observ ations dur- 
ing the 18 th and 19 th century (|Shara et"aflll985t) provide no 
indications for a second, decadal phase of high luminosity 
and low temperatures for CK Vul, as predicted by double- 
loop models, but the extinction may have hidden such an 
event. However, the low luminosity we derive is a severe 
constraint for such models. 



6.3.1 Accretion-induced pulse 

The morphology shows similarities to evolved binary sys- 
tems, such as symbiotics. This together with the compact 
radio source suggests variation on the standard VLTP sce- 
nario can be considered, involving accretion. 

A group of binary post-AGB stars, including RV Tau 
stars, show compact circumstel l ar disks, with typical s izes of 
10-100 AU (|van Winckej|2003l : Ide Ruvter et alj|2006t ). The 
disks form when a close binary (P ~ 1 yr) captures mat- 
ter lost by the central star into a circumbinary disk. The 
radio source in CK Vul has the correct size to be the rem- 
nant of such a disk. CK Vul shows similarities in particular 
with OH231.8+4.2, which shows bipolar lobes, an inner disk 
with diameter o f 50 AU as well as a 1000 AU dense torus 
|Matsuurall2006l ). 

In such a system, the post-AGB evolution may be af- 
fected by accretion from the disk on the post-AGB star 
(|Ziilstra et alj|200ll) . Such accretion could trigger a thermal 
pulse later during the post-AGB evolution than wold other- 
wise occur. The circumstellar environment can explain the 
bipolar outflows. This possibility should be explored further 
for CK Vul. 



7 CONCLUSIONS 

A compact radio source allows us for the first time to iden- 
tify the central object of the old nova CK Vul. The radio 
emission is modelled as thermal free-free emission, from a 
compact nebula of diameter ~ 60 AU. The density of the 
nebula is n c = 5 x 10 5 cm -3 . A high extinction through the 
circumstellar shell is indicated by the lack of Ha detection 
at the position of the radio source. The ionization requires 
a luminosity of the ionizing object of order lL©. 

Deep Ha images reveal a large bipolar outflow, with a 
diameter of 70 arcsec. The nebulosity is brightest near the 
position of the radio source, and shows various components. 
Comparison with previous images shows significant expan- 
sion of the nebula. Extrapolating back shows that the Ha 
structures were ejected in the AD 1670 eruption. The centre 



12 Hajduk et al. 



of expansion is, within the uncertainties, identical to the ra- 
dio source. Echelle spectra and proper motions give original 
ejecta velocities of ~ 360 km s -1 . 

The observed characteristics differ from classical novae. 
New models of nova eruption on cool, low-mass white dwarfs 
can explain some of the observational constraints. A classifi- 
cation as a Very Late Thermal Pulse appears possible, possi- 
bly accretion-induced. The radio source is proposed to be the 
remnant of a circumbinary disk, as seen in some binary post- 
AGB stars. The ejecta morphology are suggestive of binary 
interaction. The radio source provides important informa- 
tion on the environment of CK Vul. However, at present CK 
Vul still remains an enigma (|Evans et al.ll2002l ; lNa"vlor et all 
ll992l ; ISharaet al.ll 19851 ). 
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